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Abstract 

Synchronization transitions are investigated in coupled chaotic maps. Depending on 
the relative weight of linear versus nonlinear instability mechanisms associated to 
the single map two different scenarios for the transition may occur. When only two 
maps are considered we always find that the critical coupling E\ for chaotic synchro- 
nization can be predicted within a linear analysis by the vanishing of the transverse 
Lyapunov exponent At- However, major differences between transitions driven by 
linear or nonlinear mechanisms are revealed by the dynamics of the transient to- 
ward the synchronized state. As a representative example of extended systems a one 
dimensional lattice of chaotic maps with power-law coupling is considered. In this 
high dimensional model finite amplitude instabilities may have a dramatic effect 
on the transition. For strong nonlinearities an exponential divergence of the syn- 
chronization times with the chain length can be observed above ej, notwithstanding 
the transverse dynamics is stable against infinitesimal perturbations at any instant. 
Therefore, the transition takes place at a coupling e n i definitely larger than ei and 
its origin is intrinsically nonlinear. The linearly driven transitions are continuous 
and can be described in terms of mean field results for non-equilibrium phase tran- 
sitions with long range interactions. While the transitions dominated by nonlinear 
mechanisms appear to be discontinuous. 
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1 Introduction 



Synchronization is a phenomenon observed in many different contexts, rang- 
ing from epidemics spreading [1] to neuro-sciences [2]. The analysis of simple 
models pursued in the last two decades has clarified several aspects of syn- 
chronization (for a comprehensive review on the subject see Ref. [3]). 

Particularly interesting both from a theoretical and an applied point of view is 
synchronization of chaotic systems. This is a phenomenon known since many 
years [4] with important applications, e.g., in secure communications of dig- 
ital signals, laser dynamics etc. A simple but not trivial framework to study 
chaotic synchronization is represented by discrete-time dynamical systems. 
Within this context, the essence of the phenomenon can be captured already 
by considering two coupled identical maps: 



x t+1 = (l-e)f(x t )+ef(y t ) 

y t+1 =ef(x t ) + (l-e)f(y t ), (1) 

where f(x) is a chaotic map of the unit interval and e the coupling constant. 
For coupling larger than a critical value e\ a transition from a desynchronized 
state to a completely synchronized one along a chaotic orbit x* is observed, i.e. 
x l = y l for e > E\. The coupling E\ can be predicted within a linear analysis 
as the value of e for which the Transverse Lyapunov Exponent (TLE) X T , 
ruling the instabilities transversal to the line x = y, vanishes. For model (1) 
At = ln(l — 2ei) + A [3], and consequently 

1 - e_Ao 

being Ao the Lyapunov Exponent (LE) of the single (uncoupled) map. 

While in low dimensional systems the basic mechanisms ruling synchroniza- 
tion, in the absence of strong nonlinear effects, have been well established since 
long time [3,5,6], the studies devoted to synchronization in extended systems 
are still limited to few specific examples [7,8,9,10]. In the latter case, two main 
frameworks have been considered. 



On one hand, the investigation of mutual synchronization between two repli- 
cas of a spatially extended system has been pursued for systems coupled ei- 
ther via local interaction [9] or via spatio-temporal noise [8]. In both cases, 
when the amplitude of the coupling (noise) overcomes a certain threshold, syn- 
chronization is observed: after a transient, the two systems follow the same 
spatio-temporal chaotic (stochastic) orbit. 

On the other hand, another commonly observed situation corresponds to self 
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synchronization of elements belonging to the same system. This occurs for 
instance in large collections of coupled elements such as populations of neu- 
rons [11], Josephson junctions [12] or cardiac pacemaker cells [13]. In these 
systems the interaction among the elements can range from nearest neighbors 
to globally coupled. Therefore, a model able to encompass both these limiting 
cases is certainly of interest. A good candidate is represented by the following 
coupled map lattice with power-law interactions [14,15,10] 

4" = (1 - + E f{xlk) + m+k) , (3) 

which reduces to globally coupled maps (GCM's) [16] for a = 0, and to stan- 
dard coupled map lattices (CML's) with nearest neighbor coupling [17] in the 
limit a — > oo. In equation (3), t and i are the discrete temporal and spatial 
indexes, L is the lattice size (i — 1, . . . ,L), x\ the state variable, e G [0 : 1] 
measures the strength of the coupling and a its power-law decay. Since the 
sum extends up to U — (L — l)/2 the model is well defined only for odd 
L-values, and 77(a) = 2>Ylk=i k~ a is a normalization factor. Periodic boundary 
conditions are assumed. 

The synchronization between replicas of the same spatial system has been the 
subject of a more systematic study than self synchronization and is now rather 
well understood [8,18,9,19,20]. Two different mechanisms of mutual synchro- 
nization have been identified according to the predominance of linear versus 
nonlinear effects. Nonlinearly dominated transitions are observed when the 
local dynamics is ruled by a discontinuous (e.g the Bernoulli map) or "almost 
discontinuous" maps (e.g. maps possessing very high values of the first deriva- 
tives). In these linear analysis is no more sufficient to fully characterize 
the transition, because the instabilities associated with finite amplitude per- 
turbations may desynchronize the system. In both cases the transition to the 
synchronized state is typically of the second order. However, depending on the 
linear or non-linear nature of the prevailing mechanism two different universal- 
ity classes characterize the transition itself. For continuous maps (e.g. logistic 
maps) critical exponents associated with the Multiplicative Noise (MN) uni- 
versality class are usually found [21], while for (almost) discontinuous maps 
the transition belongs to the Directed Percolation (DP) class [22] . 

The aim of the present work is to clarify the effect of strong nonlinearities in the 
synchronization phenomenon, in particular in the case of self synchronization. 

As a first point, we show that strong nonlinearities may induce nontrivial ef- 
fects also in low dimensional systems such as two coupled maps. In particular, 
we find that even though the critical coupling s\ can be always predicted within 
the linear framework, the transient dynamics preceding the synchronization 
can be strongly affected by finite scale instabilities. To properly characterize 
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the latter we introduce a new indicator, the Finite Size Transverse Lyapunov 
Exponent (FSTLE), which generalizes the concept of Finite Size Lyapunov 
Exponent (FSLE) [23,24,25] to the transverse dynamics. The FSTLE extends 
the definition of TLE to finite perturbation and is therefore able to discern 
linearly from nonlinearly dominated transitions. Depending on the behavior 
of the FSTLE at finite amplitudes two different class of maps are singled out 
class I maps characterized by a decreasing FSTLE at any finite scale, and 
class II maps that present a peak in the FSTLE for some finite amplitude 
value. Moreover, in the proximity of the transition the shape of the probabil- 
ity density function (PDF), P(r), of the synchronization times r depends on 
the dynamics at finite scales and on the multifractal properties of the map. 
For class I maps two cases can be identified. For maps, such as the symmetric 
tent map and the logistic map at the crisis, which are characterized by the 
vanishing of fluctuations of the finite time LE at long times, the PDF exhibits 
a fast falloff at large r's. In particular, for the symmetrical tent map we pro- 
vide an analytical expression for P(t). For maps, such as the skew tent map, 
exhibiting modulational intermittency [3] (which is related to the persistence 
of fluctuations of the finite time LE also in the long time limit), P{t) becomes 
an inverse Gaussian distribution [26] originating by the diffusive motion of 
the (logarithm of the) perturbation in transverse space. For class II maps the 
PDF's display an exponential tail at long times which, differently from the 
previous case, is due to nonlinear effects. 



As far as spatially extended systems are concerned, Anteneodo et al. [10] have 
performed a linear stability analysis of model (3) obtaining analytically the 
critical coupling e\ for the transition. This prediction works perfectly for class 
I maps, as found for the coupled logistic maps [10]. However, as shown in 
the present paper, linear analysis may fail in class II maps. In particular, for 
large system sizes (L ^> 1) an exponential divergence of the synchronization 
times with L can be observed even for negative TLE, as a result in the limit 
L — > oo the transition takes place at a critical coupling e n \ > S\. Finite scale 
instabilities are the key elements for observing this nonlinear synchronization 
transition. Similarly to the case of diffusively coupled map lattices it is pos- 
sible to link the synchronization transition of CML with power law coupling 
with non-equilibrium phase transitions. Indeed as shown in Ref. [21] synchro- 
nization of replicas of CML with short range coupling belong either to the MN 
or to the DP universality classes. Here we shall discuss the connection of self 
synchronization of model (3) with long range spreading processes [22,27]. 



The material is organized as follows. Sect. 2 is devoted to synchronization 
of two coupled maps. In Sect. 3 coupled maps with power law coupling are 
examined. Finally in Sec. 4 we briefly summarize the reported results. 



4 



2 Synchronization of two coupled identical maps 



To analyze the synchronization of two coupled maps (1) it is useful to introduce 
the following variables: 



u 



t x l + y l 



w 



t x l — y l 



2 2 
in terms of which Eq. (1) can be rewritten as 



u 



t+1 



"'' +J ' 2 ~~ : 



f ( u t + w t ) + f{u t 



w 



(4) 



It is easily checked that the synchronized solution (corresponding to w 1 = 
and u l = f(u 1 )) is an admissible solution of (4) for any value of the coupling 
e. The stability of such a solution can be studied by considering the linear 
dynamics of an infinitesimal perturbation Sw f of the synchronized state. This 
evolves according to the linearized equation: 



5w t+1 = (1 - 2e)f'(u t )5w t 



(5) 



where /' = df/dx. Clearly the stability of the state w l = is controlled by 
the sign of the transverse (or conditional) Lyapunov exponent 



A T = ln(l - 2e) + A , 



(6) 



where Ao is the Lyapunov exponent of the single map f(x). Therefore, by 
requiring At = one obtains the critical coupling E\ (2) above which the 
synchronized state is linearly stable. Notice that for e > e h \ T coincides with 
the second Lyapunov exponent A 2 of the coupled system (1). 

The synchronization transition has been mainly studied for continuous maps, 
here denoted as class /maps, in particular for the logistic map f(x) = 4x(l— x) 
at the crisis [5] and the skew tent map [6] 



/(*) 



x/a if 0<x<« 

(1 -x)/(l -a) if a<x<l 



(7) 



It should be noticed that inside this class of maps one has to distinguish 
between two situations according to the behavior of the fluctuations of the 
finite time Lyapunov exponents in the long time limit. The vanishing of such 
fluctuations is observed for the symmetric tent map and the logistic one at the 
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Fig. 1. Typical evolution of \wt\ = \xt — yt\ starting from random initial conditions. 
From top: Bernoulli map with r = 1.5, logistic map at the crisis and the skew tent 
map with a = 2/3. For all maps, the coupling has been chosen in order to ensure 
At = — 310 -4 . The dashed straight lines display the e Ar * decay. 

crisis, while their persistence characterizes the skew tent map that represents 
the generic case [3]. 

Recent studies however pointed out that discontinuous or "almost" discontin- 
uous maps (i.e. maps with |/'| ^> 1 in some point of the definition interval), 
here termed class II maps, may give rise to non trivial interesting phenomena. 
For instance, the synchronization transition for two replicas of class II CMLs 
is not driven, as usual, by infinitesimal perturbations, but by the instability 
(spreading) of finite amplitude perturbations [8,9]. This nonlinear synchroniza- 
tion is linked to the so-called stable chaos [28]: a phenomenon characterized 
by information propagation even in the absence of chaos [29,25]. Moreover, as 
we shall show in the following also the self synchronization of a single CML 
can be strongly affected by nonlinearities. 



Representative examples of class II maps are the generalized Bernoulli map 
f(x) = rx mod 1 that is discontinuous and its continuous version 



/(*) 



a\x < x < x\ 

m(x — X2) xi < x < X2 
ai(x — x 2 ) x 2 < x < 1 



(8) 
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with a,\ = 1/xi, x\ t 2 = 1/a ± 5, m = —1/(25). As already mentioned, the 
peculiarity of such maps is that they are not only unstable with respect to 
infinitesimal perturbations but also to finite ones (for (8) this is true for suf- 
ficiently small 8 values). 

The main differences between synchronization transitions ruled by linear mech- 
anisms with respect to those driven by finite amplitude instabilities can be 
captured already by considering the evolution of The typical behavior 
of this quantity for e > S\ is reported in Fig. 1 for the Bernoulli, the logistic 
and the skew tent map. For the logistic map, apart from (short time) fluc- 
tuations due to the multiplier variations, an average linear decrease of In \w l \ 
with slope given by At is observed. For the Bernoulli map, despite the TLE 
is negative and does not fluctuate in time, exhibits instantaneous jumps 
to 0(1) values. These jumps, which are due to finite amplitude instabilities, 
are more probable for large amplitudes. Completely analogous behaviors 
have been found for the map (8), where the multipliers do fluctuate. 

As shown in Fig.l for the skew tent map with a > 1/2, \w t \ displays an 
intermittent evolution during the transient preceding the synchronization, that 
is reflected in a diffusive motion of ln|w t |. As a consequence long transients 
are observed close to the transition. This phenomenon, known as Modulational 
Intermittency [3] is induced by fluctuations of the finite time TLE that do not 
average out in the long time limit and as a consequence can be explained by 
linear analysis. 

For the sake of completeness, we mention that long synchronization transients 
have been reported also for nonlinearly coupled expanding maps [7]. In this 
model resurgences of during the transient are due to the presence of an 
invariant chaotic repelling set. We stress that for class II maps the transient 
is absolutely non chaotic i.e. their transverse dynamics in tangent space is 
contracting at any time. 

2.1 Finite Size Transverse Lyapunov Exponent 

To quantitatively characterize the different transverse-space dynamics (i.e. 
the evolution of w t shown in Fig. 1) for maps of class I and II, let us now 
introduce the Finite Size Transverse Lyapunov Exponent. The FSTLE, Ay(A), 
generalizes the concept of transverse Lyapunov exponent to finite value of the 
perturbation |u>*| = A. 

Following Aurell et al. [23] (see also Ref. [24]) we have defined the FSTLE as 
follows. We introduce a set of thresholds A n = Aok n with n — 1, . . . , N, since 
on average a transverse expanding (resp. contracting) dynamics is expected in 
the desynchronized (resp. synchronized) regime, we choose k > 1 (e.g. k = 2) 
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for e < si and k < 1 (e.g. k = 1/2) for e > E\. First, starting from a random 
initial conditions we wait for x l to relax onto its attractor. Then we initialize a 
second variable y° as y° = x° + S by choosing an initial perturbation 5 < A 
(resp. 5 > A ) if e < Si (resp. if e > £;), and let x f and y l evolve according 
to Eq. (1). Care should, of course, be taken to maintain y° within the interval 
of definition of the map. During the evolution we record the time, r(A n ,k), 
needed for \w l \ to pass for the first time from one threshold A n to the following 
one A n+1 and also the value of = w n at the moment of the passage. When 
the last threshold, A N , is reached the system is reinitialized with the above 
described procedure and this is repeated M times. The FSTLE is thus defined 
as 



where the average ([•]) is performed over the set of M different initial condi- 
tions. In the limit A — > the FSTLE converges to \ T - Notice that due to its 
definition, the FSTLE cannot measure at the same time expansion and con- 
traction rates, i.e. we have limited the analysis only to consecutive contractions 
(resp. expansions) in the synchronized (resp. desynchronized) regime. This im- 
plies that the sign of At(A) is always negative or positive in accordance with 
the investigated situation. 

In Fig. 2 A T (A) versus A is shown for different maps in the desynchronized 
regime. As expected, in all cases, for very small values of |u>*| the TLE is 
recovered. However, at larger values of |u>*| maps of the class II display an 
increase of the growth rate, i.e., A T (A) > A T (0) for some finite A value, while 
for maps of class I the FSTLE is monotonically decreasing with A. A T (A) 
provides a quantitative measure of the strength of the nonlinear effects that, 
in maps of class II, may in principle overwhelm the linear mechanisms as 
pointed out in Ref. [30,25]. 

In Fig. 3 we report the behavior of A^(A) in the synchronized regime. Again, 
while for very small A the (negative) TLE is recovered, at larger A values the 
Bernoulli-like maps display an increase which is due to the jumps 0(1) of the 
transverse perturbation. A closer inspection of the small A regime reveals that 
for the logistic map A^(A) Aj 1 only for very small scales ~ 10~ 7 , while for 
the skew tent map the asymptotic At(A) At in never exactly reached. The 
latter convergence difficulties are probably due to the intermittent behavior 
along the transverse direction, in the proximity of the transition [3]. In the 
regime e < ei the finite time effects are less important and the asymptotic 
behavior is recovered for all the investigated maps. 

The results presented in this subsection clearly indicate that finite scale insta- 
bilities can prevail on infinitesimal ones only for strongly nonlinear maps. 
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Fig. 2. At(A) vs A in the desynchronized regime for the Bernoulli map with a = 1.5 
(empty circles), its continuous version (empty boxes) with 5 = 10 , the logistic map 
(filled circles) and the skew tent map (filled triangles). The coupling constants have 
been chosen in such a way that in all systems the transverse LE is equal to At ~ 0.03 
(solid line). 
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Fig. 3. At (A) vs A in the synchronized regime for the Bernoulli map at a = 1.5 
(empty circles), the logistic map (filled boxes). The inset show the logistic map 
(filled boxes) and the skew tent map (empty boxes) . The straight lines correspond 
to the transverse LE, here At = —4 10~ 4 for all systems. 

2.2 Statistics of the synchronization times 



Notwithstanding the observed differences in the FSTLE, we found that the 
critical coupling for synchronization to occur is always given by (2) for both 
classes. This means that, at least for the case of two coupled maps, nonlinear 
effects do never modify the critical coupling. It is natural to wonder whether 
other observables related to the synchronization transition could be influenced 
by the presence of finite amplitude instabilities. As we shall show in this 
section, this is the case for the synchronization times statistics. 

Let us define the synchronization time r, as the shortest time needed to \w l \ 
for decreasing below a certain threshold 6, with the further requirement that 



9 



< for a sufficiently long successive time Tu- However, if the threshold 
value is small enough (e.g. ~ 10~ 8 — 10~ 14 ) r essentially coincides with the 
first arrival time to the considered G. The quantities of interest are the PDF's, 
P(r), and their moments (we shall focus mainly on the first moment). 



2.2.1 Class I maps 

In Fig. 4 we show P{t) for the the logistic and the symmetric tent map. The 
main feature is the presence of an exponential tail at short times and the faster 
than exponential falloff at long times. This means that synchronization times 
t ^> (r) are not observed. As a first result we derive, by following Ref. [6], an 
approximate analytical expression for P(t), which is exact for the symmetric 
tent map. 

Let us define z l = In |tu*|, and consider its linearized evolution, from Eqs. (5) 
and (6) one obtains 

z t+1 = z t + \n\f'(u t )\ + \ T -\ . (10) 

Formally the above equation can be applied only in the true linear regime, 
i.e. when — > 0, so that it is not appropriate in the early stages of the 
evolution. However, as shown in Fig. 1, for maps like the logistic one, usually 
after a reasonably short transient, the linear regime sets in and the use of 
Eq. (10) is justified for the successive evolution. From Eq. (10) it is clear that 
P(t) is related to the PDF P(z) of the variable z and to that of the local 
multipliers ln|/'(w*)|. The formal solution of Eq. (10) up to time N can be 
written as follows: 

z N — z° + X T N + A N N , (11) 

where = l/NJ2i=o,N~i l n \ f ( u% )\ ~ -V For sufficiently large N, large devi- 
ation theory tells us that the PDF of Ajv takes the formp(A) ~ exp(— Ng(A)), 
being g(A) the Cramer function [31], which is convex and has its minimum 
value at g(A = 0) = 0. It is now clear the distinction between the generic case 
in which g(A) does not collapse onto a 5-function, and the non-generic one 
in which it does, as for the symmetric tent map and the logistic one at the 
Ulam point. In the former case the dynamics of z 1 becomes a biased Brownian 
motion, with an average drift given by \ T . 

For the sake of simplicity, let us start from the tent map (Eq. (7) for a = 1/2), 
for which ln = ln(2) = A . In this case r is simply given by: 

T -ih-TTi < 12 > 

| Ay | I At | 

and P(t) can be directly related to P(z). Since, in the proximity of the tran- 
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Fig. 4. P(t) as a function of the synchronization time r : (a) tent map for 
be = (e — Ei) = 0.001, the dashed line is the analytical expression (14); (b) lo- 
gistic map for be = 0.001, the solid line refers to the expression (15) with D = 0.82 
and D' = 0.72. In all cases the PDF have been obtained by averaging over 10 6 
different initial conditions and by choosing © = 10~ 12 . 

sition, P(z = ln|iu|) assumes the form (analytically derived in Ref. [6]): 

P( z ) = J— e z e -\*F\ eZ , (13) 
| Ay | 

it can be easily obtained the related distribution 

P(r) = 2el A -l^ lne c-^ c ' ATlT+lne , (14) 



which perfectly agrees with the numerical results (Fig. 4a). 

Unfortunately the (short time) multiplier statistics of the logistic map is 
nontrivial, impeding a straightforward derivation of P(t). However, we nu- 
merically observed that P(w) ~ e~ Dw (w = exp(z)) with D almost con- 
stant in a neighborhood of the transition (namely, D = 0.82 ± 0.02 for 
-0.01 < Se < 0.01). This implies that 

P(r) = 2e D ^ +c e^ cDlXTlT+C , (15) 

which is in fairly good agreement with the numerically evaluated PDF (Fig. 4b). 
Notice that C is not simply given by D In O, and a fitting procedure is needed. 
We found C = £>' In 9, with D' = 0.72, 0.70 and 0.65 for Se = 10" 3 , 10" 4 and 
10 -5 , respectively. 

The above reported results allow us also to predict the scaling of the average 
synchronization time (r) with 5e in the proximity of the synchronization tran- 
sition. In particular, for the symmetric tent map the following expression can 
be derived 

(r>= rdrrP(r)^^hA , (16) 

J0 \At\ 
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Fig. 5. (r)) vs 5e = e — ei for the tent map (empty boxes) and the logistic one (empty- 
circles). The continuous line is the prediction (17) with 6 = 1CT 12 , the dotted one 
is obtained by a best fit of the form (r) = (a + b\n(5e))/5e, with a = 5.58 and 
b = 0.24. The data of the logistic map have been shifted by a factor 100 for plotting 
purposes. 

where [32] 

roo Q-V 00 ( — \\ n x n 

E l (x)= dy = ->y-\nx-J2 i » 

Jx y ^ nln 

being 7 ~ —0.57721. . . the Euler-Mascheroni constant. Since 20/|At| <C 1, 
approximatively we have 

Ww -j^-2 ± hwa i (17) 

where |A T | « 2e A °e = ASe (being A = In 2). Note that (17) is in perfect 
agreement with the numerical results for the tent map (see Fig. 5). For the 
logistic map, a similar dependence on 5e, namely (r) = (a + bln(Se))/Se, has 
been found. The interesting point in (17) is the logarithmic correction to the 
scaling (r) ~ Se -1 , which would have not been predicted by a naive guess 
based on (12). 

We now briefly discuss the generic situation (see Chap. 13 of Ref. [3] for 
more details). As mentioned above, z l = In performs a biased Brownian 
motion, so that computing P(r) amounts to evaluate the distribution of the 
first passage times to a threshold of a Wiener process with drift that, as a 
standard result of stochastic process theory ,is given by an inverse Gaussian 
density [26] that for long times r becomes 

P(t) ~ T- z ' 2 e~ VT with uoc\ 2 T , (18) 

where the quadratic dependence from the TLE comes from the diffusive dy- 
namics of the perturbation. 
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2.2.2 Class II maps 



For class II maps the situation is completely different. For both the continuous 
Bernoulli map (Fig. 6a) and the discontinuous Bernoulli shift (Fig. 6b) we 
observe that for e > Si the PDF's are characterized by an exponential tail at 
large r, similar to Poissonian distributions. Moreover, we also observe that the 
PDF's, once rescaled as aP(r), and reported as a function of x — (r — (r))/a 

(where a = \J (t 2 ) — (r) 2 ) collapse onto a common curve in the proximity 
of the transition. These results are particularly striking in the case of the 
Bernoulli map for which naively one would have expected a complete similarity 
with the tent map, since its multiplier f (x) = r is also constant. Here the 
presence of strong nonlinear effects makes the transient time statistics to be 
dramatically different from (14). Indeed, as seen in Fig. 1, w t is subject to 
noticeably nonlinear amplifications induced by the (almost) discontinuities in 
the map. Therefore, the expression (10) is no more appropriate to describe the 
dynamics of zt — In \ wt\, and its full nonlinear dynamics has to be taken into 
account. The latter is characterized by the fact that, with a finite probability, 
| it;* j can jump to 0(1) values at any time during the transient preceding the 
synchronization, even if the transient is not chaotic. 

This idea can be better clarified and its consequences on the P(r) better 
appreciated by considering a simple stochastic model for the dynamics of the 
transverse variable at e > e%. This model was originally proposed in Ref. [33], 
and reads as 

1 w.p. p = e AT u>* 
w t+1 =1 1 y , (19) 

I q Xt w 1 w.p. 1 — p 

where At < and w.p. is a shorthand notation for with probability. The 
underlying idea is very simple: the transverse perturbation w l is usually con- 
tracted, but with probability proportional to its amplitude can be re-expanded 
to 0(1) values, in agreement with the numerical observations. The jumps oc- 
cur when x l and y t are close but located at the opposite sides of the (almost) 
discontinuity. 

For this simple model, it is possible to derive the corresponding probability 
density P(t), which displays the same peculiar features of the PDFs reported 
in Figs. 6a and 6b. In particular, the long time tail. First, notice that the 
probability that an initial perturbation wo is never amplified to 0(1) up to 
time t is given by 



t r oo fe„-|A T |fe 

Q(t, wo) = n [1 - ^oe A H - Q(w ) = lim exp - £ °_ 



.(20) 



Interestingly, for — oo < \t < the quantity Q(wq) is always positive and 
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(t-<x>)/ct (t-<t>)/o 

Fig. 6. (a) Per as a function of (r — (t))/<7 (where a = \J (r 2 ) — (r) 2 ) for the contin- 
uous Bernoulli shift map (8) with a\ = 1.1 and 5 = 1CT 3 at different distances, Se, 
from the critical coupling. Note the fairly good collapse, indicating a common decay 
as e~ x (dashed line). The collapse is however only approximate near the peak, (b) 
Same as (a) for the Bernoulli map with r = 2 at different distances from the critical 
coupling. The dashed line gives e~ x . 

strictly smaller than 1 for any perturbation \w \ > 0, i.e. the probability that 
a perturbation of amplitude \w \ can be amplified is finite at any time. 

The PDF of the times r needed to observe w T = can be factorized as 

P e (r) = G(r-n)Q(n,l), (21) 

where G(t — n) is the probability to receive a kick at time r — n, and Q(n, 1) 
is the probability of not being amplified for the following n steps, where n is 
given by 

n = -ln6/|Ar| • 

As shown in Ref. [33] G(x) ~ exp (— ux), and therefore at large r's Eq. (21) 
can be rewritten as 

Pq(t) ocQ(n,l)exp[-i/(T + lne/|A T |)] (22) 

that confirms the Poissonian character of the PDF P(t) for the maps of class 
II. By properly normalizing Eq. (22), one obtains (r) = 1/u — Ib.Q/\\t\, 
that is in good agreement with the numerical results for the Bernoulli map (in 
particular, we considered r = 2 and = 10~ 12 ). Moreover, in the interval Se G 
[1CT 5 ; 1CT 2 ] we found that the decay rate of the PDF is directly proportional 
to the TLE, i.e v ps 0.2|A t |, this should be contrasted with the quadratic 
dependence found in the case of class I maps, see Eq. (18). 

2.3 Synchronization transition 

We conclude the investigation of two coupled maps by studying the nature 
of the transition. For this purpose let us introduce the order parameter Q(e) 
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defined as (see also Ref. [34]): 



N- 



1 

-■x N f- 



1 

•oo T 



lim — lim — > " > " \w t+Tw 



N T 

EE 

1 1=\ 



(23) 



i.e. for each value of the coupling e, one considers N different random initial 
conditions and each one is iterated for a transient T w after which the time 
average of |u>*| over a time lapse T is considered and further averaged over 
all the initial conditions. The value of the coupling giving the synchronization 
transition is then implicitly defined as fi(e n j) = 0. In principle, e n \ may dif- 
fer from Ei defined by (2), since this expression is valid only within a linear 
approximation formalism. In Fig. 7 we show the behavior of the order param- 
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Fig. 7. (a) Q(e) vs e for the Bernoulli map with r = 1.1 and its continuous version 
(8) for different values of 5. For the Bernoulli map one has that according to (2) 
£/ = 0.0454545... The order parameter is computed according to (23) averaging 
over N = 10 3 different random initial conditions, discarding the first 10 6 iteration 
and averaging over the following 10 5 iterations. O(e) has been estimated with an 
e-resolution 10 -3 . The inset shows an enlargement of the critical region done with 
a coupling resolution 10~ 4 , the same number of initial conditions and a longer 
integration transient T w = 3 10 6 iterations, (b) The same for the logistic map at 
the crisis. Here (2) predicts e\ = 0.25. The blow up of the critical region (in the 
inset) shows that the transition is continuous. Note also that for 8e <C 1 we observe 
W oc 6e. 

eter (23) as a function of e for both classes of systems. Two observations are 
in order. First, the transition to the synchronized state always occurs at the 
critical coupling defined by (2) (i.e. e n \ = s\). Second, the transition is con- 
tinuous for the Logistic map (simulations show that also the tent map, which 
belongs to the same class, has a continuous transition) and discontinuous for 
the Bernoulli shift. This confirms the results reported in Ref. [34]. However, 
since the transition for the continuous version of the Bernoulli shift map (8) 
is steep but clearly continuous (see the inset of Fig. 7a), this does not seem to 
be a generic property for all the maps of class II. 



The behavior of the Bernoulli shift map can be easily understood by consid- 
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ering the map controlling the dynamics of w t , i.e. 



w 



t+i 



[1 - 2e)rw t 



x l e h y* e h 
(1 - 2e) (rw* + 1) x* e I , y* e h 
(1 - 2s) (rw l - 1) x l e h , y l e h 



(24) 



being Iq = [0 : 1/r] and I\ = ]l/r : 1]. From the above expression, it is 
evident that the attractor in the plane has always a finite width 

along the transversal direction, unless x* — y*. This explains the discontinuity 
observed at S\. In the continuous map (8) the attractor in the plane (w*, w t+1 ) 
has a "transverse" width that decreases continuously to zero for e — > E\. In 
conclusion, the discontinuity in the transition is a pathology of the Bernoulli 
map, and possibly of other discontinuous maps. 

The order parameter (23) can be seen as the time average, once discarded an 
initial transient, of the following quantity 



wit) = (K|> 



fl(e) for t — > oo) , 



(25) 



where < • > indicates the average over many different initial conditions. Once 
w is initialized 0(1), depending on whether e is smaller or larger than Ei two 
different asymptotic behavior are observed. In the desynchronized regime, ob- 
viously Wit) goes to a finite value Q(e), while above the synchronization tran- 
sition Wit) ^0 with a decay law determined by the nature of the considered 
maps. For class I maps, the decay is ruled by the TLE: W(t) ~ exp(A^t) with 
At < 0. For class II maps W(t) shows an initial exponential decay ~ exp(— /it) 
(with /i > 0), due to the effect of the resurgences averaged over many different 
initial conditions, followed by a final linear decay exp(\ T t). 



The nonlinear rate fi can be simply related to the exponential tail of the PDF 
of the first arrival times, by assuming (as indeed observed) the following decay 
Pq(t) oc exp [—v(t — jf|j)]- Since W(t) is the average amplitude value of the 
w l at time t, one has 



W(t)= CdQ P e (t) oc e~ vt /"dG 6 e~"TX% = e- vt CdQ 1 ^ . (26) 
Jo Jo Jo 

The above result tells us that /i — v (if the integral does not diverge, i.e if 
v < | At- I ) ■ This is confirmed by numerical checks. 
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3 Coupled maps with power law coupling 



In this section we analyze high dimensional systems, namely the power-law 
coupled maps defined in Eq. (3). In particular, for class II maps we shall 
show that there exist situations where, due to the nonlinearities, the synchro- 
nization time diverges exponentially with the number of maps. As a conse- 
quence, the transition takes place at a (nonlinear) critical coupling e n i larger 
than the linear value Si. A similar phenomenon has been observed for GCM's 
with nonlinear coupling, where the synchronization time divergence is due to 
a chaotic transient [7]. Chaotic transients, diverging exponentially with the 
number of coupled elements, have been reported also for spatially extended 
reaction-diffusion systems [35] and for diluted networks of spiking neurons [36]. 
However, the emphasis of our work is on non-chaotic transients similar to the 
stable-chaos phenomenon [28]. 

Let us start by reviewing the linear theory developed in Ref. [10], which is 
able to account for the synchronization transition of class I maps. 



3.1 Lyapunov Analysis 



When nonlinear effects are not sufficiently strong, excluding the pathological 
cases of chaotic transients, the critical coupling for observing the synchro- 
nization transition can be predicted by computing the Lyapunov spectrum. 
In particular, since above the synchronization transition the transverse Lya- 
punov exponent coincides with the second Lyapunov exponent, A 2 , it suffices 
to evaluate the dependence of the latter on e. 

In order to compute the Lyapunov spectrum of the model (3), it is necessary 
to consider the tangent space evolution: 



5x\ +1 = (1 - e)f'(xl)5xl + V f'W-k)te\-k + f\x\ +k )5x\ +k 



In Ref. [10] it has been shown that the full Lyapunov spectrum can be easily 
obtained for the Bernoulli map: 



Afc = lnr + In 



1 - e + —r-rh 
77(a) 



(28) 



where 



^^c^H-lWL) t = l,...,L. 



(29) 
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In this notation the maximal Lyapunov exponent is given by Ai = lnr and 
the second Lyapunov exponent is A 2 . 

For generic maps it is still possible to obtain the Lyapunov spectrum in an 
analytical way, but only in the synchronized state, by substituting lnr with 
the maximal LE of the single uncoupled map A in (28). 

From a linear analysis point of view one expects that the synchronization 
transition has to occur when A 2 = 0. Therefore, the following expression for 
the critical line in the (a,e)-plane can be derived: 

As shown in Ref. [10] this prediction is well verified for the logistic map, which 
is here used as a benchmark for our codes. Let us conclude this short review 
by noticing that in the limit L — > oo synchronization can be achieved only for 
a < 1 [10]. Therefore we shall limit our analysis to this interval. 

3.2 Nonlinear synchronization transition 

Let us now study the critical line for class II maps, exemplified by the Bernoulli 
shift map and its continuous version. First of all we introduce the main ob- 
servables and the numerical method employed to determine it. 

A meaningful order parameter for the transition is represented by time average 
of the following mean field quantity : 

S(i) = yX><-^|, ^ = yXX (31) 
L i=i L i=i 

This can be operatively defined as follows : firstly the system (3) is randomly 
initialized and iterated for a transient time T w proportional to the system size 
L, then the time average of S(t) is computed over a time window T, i.e. (S)t = 
l/TJ2t=i,T Finally the state of the system is defined as synchronized if 
(S)t < ©, being a sufficiently small value (10~ 8 — 10~ 10 is usually enough). 
The coupling value corresponding to the synchronization transition is then 
obtained by using a bisection method: chosen two coupling values across the 
transition line one corresponding to a desynchronized state (e^) and the other 
to a synchronized case (e 8 ) a third value is selected as e m = (e d + e s )/2. If 
at this new coupling value the system synchronizes (resp. not synchronizes) 
e rn is identified with the new e s (resp. e^). The procedure is then repeated 
until (e s — €d) < e (in our simulations 9 e = 10 -3 — 10~ 5 ). Finally the critical 
coupling is defined as e n i = (e s + e<f)/2. The algorithm, tested on the logistic 
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map, was able to recover (30) with the required accuracy. In general one has 
that, once fixed a, e n \ will be a function of T w for a given L. 



0.4 
0.3 

w 

0.2 
0.1 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

a 

Fig. 8. Synchronization transition lines for the Bernoulli shift map with r = 1.1 at 
L = 17, 33, 65, 129, 257. The solid lines are the analytically estimated "linear" 
transition values Eq. (30), while the symbols refer to the numerically obtained 
values. We used T w = 10 5 L,T= 10 3 L, 6 = 10~ 8 and 6 e = 10~ 3 . 

In Refs. [7,10] it has been employed as a synchronization indicator the time 
average of the following mean-field quantity: 

(32) 

We have verified that this indicator gives results completely analogous to (S) T 
in all the considered cases. 

In Ref. [10] the authors have reported for logistic coupled maps a very good 
agreement between q given by (30) and e n \ estimated via (-R')t- However, 
this is not the case for the Bernoulli map, as shown in Fig. 8. In this case 
(depending on the slope of the map r, on a and on the chain length L) 
strong disagreements between the linear transition line given by e\ and the 
numerically obtained values are observed. These disagreements are typical of 
class II maps. Since the nonlinear effects locally desynchronize the system even 
if A 2 < 0, in general e n i > e x . 

In obtaining the data reported in Fig. 8, the cpu time restrictions forced us 
to employ a large but somehow limited transient time (namely, T w = 10 5 L). 
Therefore, we checked for the dependence of the results on T w for fixed L. In 
particular, we measured e n i for two a-values only (namely, a = 0.3 and 0.8) for 
several chain lengths and transient times. This analysis has been performed for 
the logistic map at the crisis, for the Bernoulli shift map with r — 1.1 and for 
its continuous version. The results are reported in Fig. 9. For the logistic map 
we found that, for any value of L and a, a relatively small value of T w (xs 10 4 ) 
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Fig. 9. e n l as a function of L for different T w in 10 3 — 10 7 for the continuous Bernoulli 
shift map with r = 1.1 and 5 = 3 x 1CP 4 at a = 0.3 (a) and a = 0.8 (b). The same for 
the Bernoulli shift with r = 1.1 for a = 0.3 (c) and for a = 0.8 (d). The continuous 
solid lines are the theoretical values obtained in the linear analysis framework (30). 

was sufficient to observe a clear convergence of e n i to the linear value e\. This 
is also the case of the Bernoulli map and its continuous version for a = 0.3 
(though as one can see in Fig. 9c the discontinuous map is characterized by 
a slightly slower convergence than its continuous version). On the contrary 
for a = 0.8 with L > 21 we were unable to observe synchronization even for 
T w = 10 7 L, while at T w ps 10 4 L the logistic map had already converged to e\. 
Summarizing for r = 1.1 we found that for a < 0.5 the Bernoulli shift map and 
its continuous version always converge toward the linear critical value, while 
for a > 0.5 at increasing L the critical coupling e n i becomes more and more 
independent of the transient time T w (see the right of Fig. 9 for L = 129). 
However, for a smaller value of r and for large L (we investigated the case 
r = 1.01 and L = 101,501) we observed e n i > E\ for any value of a G [0; 1]. 
In particular, we stress that also in the globally coupled case (corresponding 
to a = 0) we have clear evidences that e n i > e t . This is due to the fact that 
for r — > 1 nonlinear finite amplitude instabilities becomes more and more 
predominant with respect to the mechanism of linear stabilization (see [30,29] 
for more details). 

On the basis of the previous results, it is natural to conjecture that in the 
limit L — > oo the non-linear transition is well denned, i.e. that the limit 



e„*(a) = l im r lim e n i(a,L,T w ) 



(33) 



exists and is typically larger than ej(a) = Hindoo ej(a, L). Note that in the 
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above expression the order of the two limits is crucial, we expect that perform- 
ing at fixed L the limit T w — > oo we should always observe a convergence to e\ . 
However, as shown in the following, the times to reach synchronization may 
diverge exponentially fast with L in the region just above E\ making rapidly 
infeasible this limit. 

3.3 Synchronization Times 

As for the case of two coupled maps, we study now the synchronization time 
statistics. In particular, we consider the system above the linear transition line 
e > Si and measure the first passage times r needed for S(t) (32) to decrease 
below a given threshold 6. In this way we determine the corresponding PDF 
and the related momenta. 

In the high dimensional case, it is fundamental to analyze the dependence 
of the synchronization time on the the system size L. However, as reported 
in Eq. (30), the critical value Si itself depends on L. Therefore to perform a 
meaningful comparison of systems of different sizes we considered situations 
characterized by the same linear behavior in the transverse space, i.e. having 
the same value of A 2 . This means that for each length L we chose the coupling 
strength according to 



Further tests performed at different lengths but at a fixed distance As = 
e — Si(L) from the critical line give essentially the same results. 

As shown in Fig. 10a, for the logistic map (but these results can be extended to 
all the maps belonging to class I, see Fig. 10b) the PDF's of the synchronization 
times obtained at constant A2 display a very weak (almost absent) dependence 
on the system size, and are qualitatively similar to that found for two coupled 
maps (compare Fig. 10a with Fig. 4). Moreover, measurements of the average 
synchronization times (r), done by fixing the "distance" Ae from the critical 
line, exhibit a clear tendency to saturate for increasing L (Fig. 11). Therefore, 
in the limit L — > 00 the synchronization time will not diverge. 

For the maps of class II the situation is different. Here, as shown in Fig. 10c, 
for values of a sufficiently small P(r) is weakly dependent on the system size 
as found for the logistic map. On the other hand, by considering more local 
couplings (e.g. a = 0.8 in Fig. lOd), the tail of the PDF becomes more and 
more pronounced as the system size increases. Results obtained by fixing the 
distance from the critical coupling instead of the value of the TLE display 
qualitatively similar features. 




(34) 
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Fig. 10. -P(t) vs r for different system sizes (see label) at fixed A2 = —0.01 for 
different maps: (a) logistic map at the crisis with a = 0.8 (for a = 0.3 we obtained 
qualitatively similar results); skew tent map (7) with a = 2/3 (b); Bernoulli shift 
map with r = 1.1 for a = 0.3 (c) and a = 0.8 (d). In the latter case sizes larger 
than L = 17 were not drawn for the sake of clarity of the plot. The PDF's have 
been obtained in all cases by considering 10 4 different initial conditions and with 

e = kt 12 . 
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Fig. 11. Average synchronization time, (r), as a function of the system size L for 
the logistic map at the crisis with a = 0.6 estimated at a fixed distance from the 
linear threshold Ae = 0.05. Data have been averaged over 10 5 initial conditions 
with G = 10~ 12 . 

This picture is further confirmed by examining (r) as a function of L for 
fixed A2 (i.e. by choosing e according to (34)). As one can see in Fig. 12, (r) 
displays a dramatic dependence on L. In particular, at A2 > —0.01 (see the 
inset of Fig. 12) (r) grows exponentially with L, while a power-like scaling is 
observable for A 2 ~ —0.05, at least for the chain lengths we could reach. Deeper 
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L 

Fig. 12. Average synchronization times (r) as a function of L for the Bernoulli map 
with r = 1.1 and a = 0.8 for various values of the second Lyapunov exponent (see 
label). The inset shows the case of A2 = —0.01 in linear scale. Data have been 
averaged over 10 4 initial conditions and = 10~ 12 . 

inside the "linear" synchronization region, i.e. for A2 = —0.1, we observed a 
saturation of (r) with the system size. It is worth stressing here that for 
a = 0.3 and r = 1.1, where no appreciable discrepancies between e n i and Ei 
have been observed, (r) saturates for large L. Nevertheless the corresponding 
PDF exhibits the usual tail at long times characteristic of maps of class II but 
with weak dependence on L (see Fig. 10c). Notice also that, in principle, the 
transition from the exponential dependence to the saturation can be used to 
estimate e n i (see Ref. [8]), however for the present model a systematic study 
of this aspect is infeasible due to the required computational resources 

We thus found clear indications that a nonlinear synchronization transition 
can be observed for class II maps when finite size nonlinear instabilities are 
sufficiently strong to overcome the effects of stabilization associated with the 
local linear dynamics and to the spatial coupling. This suggests that nonlinear 
effects tend to decouple the single units of the chain. This decoupling can be 
indeed related, by means of the following simple argument, to the observed 
exponential divergence with L of the synchronization times at e n \ > e > E\. 
Synchronization happens when the difference s\ = \x\ — x t \ is contracted for 
ever in each site until S(t) = s\ decreases below threshold. The probability 
that for a single site the initial difference s° = wo will be never amplified to 
0(1) is given by Q(wq) < 1 (see Eq. (20)), therefore by assuming that each 
sites is completely decoupled from its neighbors the probability of contraction 
for the whole chain is given by p = Q L {1) = exp [LlnQ(l)] (for simplicity we 
set wo = 1 without loss of generality). And the probability that the system 
settles onto this "contracting" state after n steps is 

P n ~ (1 — p) n ~ 1 p ~ exp(— np) (35) 

this quantity represents a good approximation of P{j). Therefore one expects a 
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Poissonian decay for the PDF of the synchronization times with an associated 
average time < r >= 1/p oc exp [Lln(l/Q(l))] exponentially diverging with 
the chain length, as indeed observed. 



3.4 Properties of the transitions 

Let us now characterize the synchronization transition in the spatially ex- 
tended model (3) within the framework of non-equilibrium phase-transitions. 
A similar parallel was recently established in a series of works concerning syn- 
chronization of two replicas of CML's with nearest neighbor coupling [8,9,33,37] 
In these studies it has been shown that the synchronization transition is con- 
tinuous and belongs to the Multiplicative Noise (resp. Directed Percolation) 
universality class depending on the linear (resp. nonlinear) nature of the pre- 
vailing mechanisms. In both cases one observes a transition from an active 
phase (characterized by (S)t = So > 0) to a unique absorbing state (identi- 
fied by (S)t = 0). We indicated with S the order parameter, however for sake 
of clarity it should be said that for two replicas of a CML this corresponds to 
S(t) = J2i=i,L \ x \ —y\\ an d is not given by expression (32) used in the present 
paper for characterizing self synchronization. The rationale for using the same 
symbol is that the two definitions embody essentially the same information: 
< S >t is a measure of the density of non synchronized (active) sites. 

In the proximity of the transition (if continuous) the scaling behavior of the 
saturated density So of active sites is given by 



while at the critical point e = e c the density of active sites scales as S(t) ~ t 
and the average synchronization time diverges as 



We denoted as e c the critical coupling to avoid at this stage distinctions be- 
tween linearly and nonlinearly driven synchronization. Since continuous non- 
equilibrium phase transitions are typically characterized by three independent 
critical exponents, once (5, [3, z) are known all the other scaling exponents can 
be derived [22]. 

Let us now analyze CML's with power-law coupling (3). We start with class 
I maps, for which the transition is completely characterized by the linear 
dynamics, i.e. e c = E\. In particular, we consider the logistic maps at the Ulam 
point for three different a-values (namely, a = 0.1,0.3 and 0.6). The first 
observation is that the critical properties of the model seem to be independent 
of a: we measured exactly the same critical exponents for the all a values, 



So ~ (e c — eY for e < e, 



(36) 



< t >~ L z . 



(37) 
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Fig. 13. Temporal evolution of the order parameter S(t) for e = e c — 10 -3 (filled 
circles) e = e c (empty boxes) and e = e c + 1CT 3 (crosses). The solid straight line 
displays the power law decay Data refer to the logistic map with a = 0.1 and 
L = 501 and are obtained after averaging over 1, 500 — 3,000 initial conditions. In 
the inset the average value of the order parameter (S)t as a function of e c — e is 
reported. The straight line shows the scaling behavior (e c — e) 1 . Data refer to the 
logistic map with a = 0.3 and L = 501 the average is performed over 1, 000 different 
initial conditions. 

namely S ~ 1, f3 ~ 1 (see Fig. 13) and z ~ (see Fig. 11). These exponents 
coincide with the mean-field exponents reported in Ref. [38] for a model of 
anomalous Directed Percolation (DP). The model differently from standard 
DP considers the spreading of epidemics in the case of long-range infection, 
namely in one spatial dimension the probability distribution for a site to be 
infected at distance r decays as l/r 1+cr . The critical behavior of this model 
can be obtained by considering a Langevin equation with power-law decaying 
spatial coupling for the coarse grained density of infected sites. In Ref. [38], the 
authors found that the critical exponents vary continuously with a, but below 
a critical value (a c = 1/2 in Id) the exponents coincide with the corresponding 
mean-field results, namely (3mf — 8mf — 1 and zmf = In the limit a — > 
the latter exponents are identical to the ones we have found for the logistic 
coupled maps. In the following we shall give an argument to explain these 
similarities and differences. 



As shown in Refs. [8,9,33,37] there is a deep connection between the synchro- 
nization problem for two replicas of a diffusively coupled CML and nonequi- 
librium phase transitions. Indeed there the difference field d\ = \x\ — y\\ can be 
mapped onto the density of infected (active) sites and an appropriate Langevin 
equation describing the evolution of d\ can be derived in proximity of the tran- 
sition [20]. In our case of CMLs with power-law coupling we expect that the 
corresponding Langevin equation for the spatio-temporal coarse grained de- 
fect density s\ = \x\ — x t \ should contain a long-range interaction with a power 
law spatial coupling decaying with an exponent a. Therefore, in the proximity 
of the transition, it seems natural to map the dynamics of (3) onto the model 
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Fig. 14. Temporal evolution of the averaged order parameter S(t) for the Bernoulli 
map with r = 1.1, a = 0.8, L = 501, and different coupling strength above and 
below e n i (see the labels). Average is done over 100 — 500 initial conditions, tests 
on shorter times and more initial conditions confirm the reliability of the statistics. 

studied in [38] once the exponent a is identified with 1 — a. However, in (3) 
the coupling is rescaled by the factor i](a) ~ L 1 ^ a to avoid divergences in the 
limit L — > oo. The rescaling amounts to consider effective interaction 

on large scales of the type 1/r independently of a. This may explain the fact 
that the exponents characterizing the transition of model (3) coincide with 
the mean- field values reported in [38] for a = 0. A similar rescaling was per- 
formed in Ref. [39] to show, for a chain of power-law coupled rotators, that 
all the equilibrium properties of the system coincide for < a < 1, once suit- 
ably scaled. Let us also remark that since we are dealing with a mean-field 
case the nature of the noise entering in the Langevin equation, which distin- 
guishes Directed Percolation from Multiplicative Noise, is irrelevant [21,22]. 
This leads us to conclude that, for continuous synchronization transitions, we 
would not expect any differences in the measured exponents between linearly 
and nonlinearly driven transitions. 

In examining the critical properties associated to coupled class II maps, we 
confronted with hard numerical difficulties related to the fact that e c = e n i 
cannot be measured with the required precision. Moreover, due to the presence 
of exponentially long transients in the proximity of the nonlinear transition, 
simulations of model (3) with class II maps become extremely time consuming. 
For these reasons we were unable to find conclusive results. However, within 
these limitations, our analysis (performed up to L = 501 and up to integration 
times ~ 10 6 ) suggests that for class II maps (both the Bernoulli map and its 
continuous version (8) with r = 1.1 were considered) the transition may be 
discontinuous as suggested from the data shown in Fig. 14. In the figure the 
temporal evolution of the order parameter S(t) averaged over many initial 
conditions is reported in proximity of e n \ and no indication of critical behaviour 
is observable. This indicates that in class II CMLs: short range interactions 
give rise to DP-like continuous transitions, while power-law coupled maps 
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exhibit discontinuous transitions. We should mention that the introduction 
of long-range interactions in DP model can similarly modify the nature of 
the transition from continuous to discontinuous, as shown in a recent paper 
[27]. The authors have introduced a generalized directed percolation model, 
in which the activation rate of a site at the border of an inactive island of 
length t is given by 1 + a/£ a . In particular, in Ref. [27] it has been shown that 
the transition is continuous for a > 1 with exponent coinciding with DP, and 
discontinuous for < a < 1. 



4 Conclusions 



In the present paper we have examined the influence of strong nonlinearities in 
the transverse synchronization transition of low and high dimensional chaotic 
systems. In particular, for two coupled maps we have shown that finite am- 
plitude nonlinear instabilities can give rise to long transients, preceding the 
synchronization, even when the dynamics is transversally stable at any in- 
stant. This should be contrasted with the results found for continuous maps 
where long transients may happen only as a result of transient chaotic or in- 
termittent transverse dynamics. In high dimensional systems with long range 
(power-law) interactions strong nonlinearities may invalidate the linear crite- 
rion to locate the critical coupling. In this case the transition occurs, due to 
nonlinear mechanisms, at a larger coupling value. The nonlinear transition is 
characterized by the emergence of transients diverging exponentially with the 
system size even above the linear critical coupling. The origin of this transition 
is closely related to the stable chaos phenomenon. The synchronization phe- 
nomena both in linearly and nonlinearly driven systems have been compared 
with models for long ranged contact processes. We found that, for linearly 
driven systems, the transition is continuous and the critical exponents are 
given by a mean field prediction. For nonlinearly driven systems, though the 
results are not conclusive, evidence of a discontinuous transition have been 
found. 
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